Quantum mechanics and low energy nucleon dynamics 
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We discuss the problem of consistency of quantum mechanics as applied to low energy nucleon 
dynamics with the symmetries of QCD. It is shown that the dynamics consistent with these symme- 
tries is not governed by the Schrodinger equation. We present a new way to formulate the effective 
theory of nuclear forces as an inevitable consequence of the basic principles of quantum mechanics 
and the symmetries of strong interactions. We show that being formulated in this way the effective 
theory of nuclear forces can be put on the same firm theoretical grounds as the quantum mechanics 
of atomic phenomena. In this case the effective theory allows one to describe with a given accuracy 
not only two-nucleon scattering, but also the evolution of nucleon systems, and places the constraints 
on the off-shell behavior of the two-nucleon interaction. In this way we predict the off-shell behavior 
of the S wave two-nucleon T-matrix at very low energies when the pionless theory is applicable. 
Further extensions and applications of this approach are discussed. 



I. INTRODUCTION 

Ideas from the foundations of quantum mechanics are 
being applied now to many branches of physics. In the 
quantum mechanics of particles interacting through the 
Coulomb potential one deals with a well-defined inter- 
action Hamiltonian and the Schrodinger equation gov- 
erning the dynamics of the theory. This theory is per- 
fectly consistent and provides an excellent description of 
atomic phenomena at low energies. It is natural to ex- 
pect that low energy nuclear physics can be described 
in the same way. However, one has not yet constructed 
a fundamental nucleon-nucleon (NN) potential. Nowa- 
days there exist phenomenological NN potentials which 
successfully describe scattering data to high precision, 
but they do not emerge from QCD and contain ad hoc 
form factors. A first attempt to systematically solve the 
problem of low energy nucleon dynamics and construct 
a bridge to QCD was made by Weinberg 0. He sug- 
gested to derive a NN potential in time-ordered chiral 
perturbation theory (ChPT). However, such a potential 
is singular and the Schrodinger (Lippmann-Schwinger) 
equation does not make sense without regularization and 
renormalization. This means that in the effective field 
theory (EFT) of nuclear forces, which following the pio- 
neering work of Weinberg has become very popular in nu- 
clear physics (for a review, see Ref. |2j), the Schrodinger 
equation is not valid. On the other hand, the whole for- 
malism of fields and particles can be considered as an 
inevitable consequence of quantum mechanics, Lorentz 
invariance, and the cluster decomposition principle 
Thus in the nonrelativistic limit QCD must produce low 
energy nucleon physics consistent with the basic princi- 
ples of quantum mechanics. However, as it follows from 
the Weinberg analysis, QCD leads through ChPT to the 
low energy theory in which the Schrodinger equation is 
not valid. This means that either there is something 
wrong with QCD and ChPT or the Schrodinger equa- 
tion is not the basic dynamical equation of quantum the- 
ory. Meanwhile, in Ref. Q it has been shown that the 
Schrodinger equation is not the most general equation 



consistent with the current concepts of quantum physics, 
and a more general equation of motion has been derived 
as a consequence of the basic postulates of the Feynman 
and canonical approaches to quantum theory. Be- 
ing equivalent to the Schrodinger equation in the case 
of instantaneous interactions, this generalized dynamical 
equation permits the generalization to the case where the 
dynamics of a system is generated by a nonlocal-in-time 
interaction. The generalized quantum dynamics (GQD) 
developed in this way has proved an useful tool for solv- 
ing various problems in quantum theory Q . 

The formalism of the GQD allows one to consider the 
problem of consistency of quantum mechanics with the 
low energy predictions of QCD from a new point of view. 
From this viewpoint, in investigating the consequences 
of ChPT we must not restrict ourselves to the assump- 
tion that the NN interaction can be parameterized by a 
NN potential, and low energy nucleon dynamics is gov- 
erned by the Schrodinger equation. This dynamics may 
be governed by the generalized dynamical equation with 
a nonlocal-in-time interaction operator when this equa- 
tion is not equivalent to the Schrodinger equation, and 
hence the above divergence problems may be the cost of 
trying to describe low energy nucleon dynamics in terms 
of Hamiltonian formalism while this dynamics is really 
non-Hamiltonian. In the present paper, by using the ex- 
ample of the ^Sq channel, we will show that from the 
analysis of time-ordered diagrams for the T-matrix in 
ChPT it really follows that nucleon dynamics at very low 
energies is governed by the generalized dynamical equa- 
tion with a nonlocal-in-time interaction operator. The 
GQD provides a new way to formulate the effective the- 
ory of nuclear forces as an inevitable consequence of the 
basic principles of quantum mechanics and the symme- 
tries of QCD, and, as it will be shown, in this way one 
can construct the theory that is perfectly consistent and 
free from ultraviolet (UV) divergences. 

The aim of the present paper is twofold. From the 
one hand, we intend to show that the low energy pre- 
dictions of QCD are consistent with the basic principles 
of quantum mechanics, but in this case a new insight 
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into these principles provided by the formahsm of the 
GQD is needed. On the other hand, it is our intention to 
demonstrate that this insight into the basic principles of 
quantum mechanics opens new possibilities for describ- 
ing low energy nucleon dynamics. In Sec. II we briefly 
consider the main features of the formalism of the GQD 
developed by one of the authors (R.G.) in Ref. ^. We 
mainly focus on the physical meaning of the generalized 
dynamical equation which is a direct consequence of the 
principle of the superposition of the probability ampli- 
tudes and the requirement that the evolution operator 
be unitary. This equation plays a key role in the present 
work. 

In Sec. Ill we consider the pionless theory of nuclear 
forces. By focusing the attention on the ^Sq channel, 
we show that the requirement that the two-nucleon {2N) 
T-matrix satisfy the generalized dynamical equation and 
be consistent with the symmetries of QCD completely de- 
termines its form. This means that the pionless effective 
theory can be formulated as an inevitable consequence 
of the symmetries of QCD and the basic principles of 
quantum mechanics, if this principles are entered into the 
theory via the generalized dynamical equation. In Sec. IV 
we show that, starting from the expression for the 2N T- 
matrix obtained in this way one can organize calculations 
of observables in the spirit of the effective theory, i.e., by 
expanding amplitudes in powers of Q/A where Q being 
some low energy scale, and A is the scale at which the 
theory is expected to break down. This makes it possible 
to calculate amplitudes up to some order in Q/A keep- 
ing the terms in this expansion up to the corresponding 
order. In this way we reproduce all results for the NN 
scattering amplitudes obtained in the standard EFT of 
nuclear forces. It is important that these results are re- 
produced starting with a well defined 27V T-matrix with- 
out resorting to the regularization and renormalization 
procedures. This shows that really the ChPT does not 
lead, in the nonrelativistic limit, to a low energy theory 
with UV divergences, i.e., low energy predictions of QCD 
are consistent with the principles of quantum mechanics. 

In Sec.V we show that the use of the generalized dy- 
namical equation allows one to formulate the effective 
theory as a perfectly consistent theory free from UV di- 
vergences. Being formulated in this way, the effective the- 
ory keeps all advantages of the traditional nuclear physics 
approach. At the same time, its advantage over the tradi- 
tional approach is that it allows one to find constraints on 
the off-shell behavior of the 2N T-matrix that are placed 
by the symmetries of QCD. In contrast, the realistic po- 
tentials that has been developed within the traditional 
approach cannot guarantee a reliable off-shell extrapola- 
tion of the 2N T-matrix since they are all constrained by 
the 2N phase shift analysis. 

In Sec. VI we show that, in despite of the wide-spread 
opinion, from the first principles of quantum mechanics it 
follows that in general renormalization is not necessary to 
separate the low energy physics that is described by the 
effective theory of nuclear forces from the underlying high 



energy physics. An important feature of the generalized 
dynamical equation is that it provides such a separation 
without renormalization. The only thing that is needed 
for this is that at low energies nucleons emerge as the 
only effective degrees of freedom. 

In Sec. VII the proposed formalism is considered from 
the point of view of the Weinberg program for physics of 
the two-nucleon systems. We show that our formulation 
of the effective theory of nuclear forces can be considered 
as a new way of the realization of the Weinberg program. 
This way leads to a perfectly consistent effective theory 
that allows one not only to calculate the 27V scattering 
amplitudes but 27V T-matrix and hence the evolution and 
Green operators. Moreover, the effective theory devel- 
oped in this way allows one to determine the off-shell be- 
havior of the 27V T-matrix as an inevitable consequence 
of the first principles of quantum mechanics and the sym- 
metries of QCD. This is very important, because, as is 
well known, the off-shell behavior of the 27V T-matrix 
may play a crucial role in solving the many-nucleon prob- 
lem and is an important factor in calculating in-medium 
observables 8] and in microscopic nuclear structure cal- 
culations. This results, for example, in the fact that the 
predictions by the Bonn potential for nuclear structure 
problems differ in a characteristic way from the ones ob- 
tained with local realistic potentials Q . The off-shell am- 
biguities of realistic 7V7V potentials are argued to be one 
of the main causes of many problems in describing three- 
nuclcon systems. One of such problems is the Ay puzzle 
which refers to the inability to explain the nucleon vec- 
tor analyzing power Ay in elastic nucleon-deuteron (Nd) 
scattering at low energy. As has been shown in Ref. 0| , 
it is not possible with reasonable changes in realistic 7V7V 
potentials to increase the Nd Ay and at the same time 
to keep the 27V observables unchanged. The same situa- 
tion also takes place in the case of chiral potentials [lTI |. 
This means that the introduction of three-nucleon forces 
(3NF) is needed for resolving the problem. However, as 
it turned out, conventional 37VF's change the predictions 
for Nd Ay only slightly and do not improve them [12- 
14]. This motivated new types of 37VF's [15-17]. On the 
other hand, reliable 37V calculations and even testing for 
the presence of 37V forces require to constrain the off- 
shell properties of the 27V. Moreover, these properties 
play a crucial role in a new 37VF proposed by Canton 
and Schadow |0. These aspects of the three-nucleon 
problem are discussed in Sec. VI. We show that the fact 
that the GQD allows one to predict the off-shell behavior 
of the 27V T-matrix, by using the symmetry constrains 
placed by a QCD, may open new possibilities for solving 
three-body problem in nuclear physics. 



II. GENERALIZED QUANTUM DYNAMICS 

The basic concept of the canonical formalism of quan- 
tum theory is that it can be formulated in terms of 
vectors of a Hilbert space and operators acting on this 
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space. In this formalism the postulates that establish 
the connection between the vectors and operators and 
states of a quantum system and observables are used to- 
gether with the dynamical postulate according to which 
the time evolution of a quantum system is governed by 
the Schrodinger equation. In the Feynman formalism 
quantum theory is formulated in terms of probability 
amplitudes without resorting to vectors and operators 
acting on a Hilbert space. Feynman's theory starts with 
an analysis of the phenomenon of quantum interference. 
The results of this analysis which leads directly to the 
concept of the superposition of probability amplitudes 
are summarized by the following postulate i5|: 

The probability of an event is the absolute square of a 
complex number called the probability amplitude. The 
joint probability amplitude of a time-ordered sequence of 
events is the product of separate probability amplitudes 
of each of these events. The probability amplitude of an 
event which can happen in several different ways is a sum 
of the probability amplitudes for each of these ways. 

According to this assumption, the probability ampli- 
tude of an event which can happen in several different 
ways is a sum of contributions from each alternative way. 
The Feynman formulation is based on the assumption 
that the history of a quantum system can be represented 
by some path in space-time, and hence the probability 
amplitude of any event is a sum of the probability ampli- 
tudes that a particle has a completely specified path in 
space-time. The contribution from a single path is postu- 
lated to be an exponential whose (imaginary) phase is the 
classical action (in units of h) for the path in question. 
This assumption is not as fundamental as the above prin- 
ciple of the superposition of probability amplitudes which 
follows directly from the analyzes of the phenomenon of 
quantum interference. This fact is emphasized, for ex- 
ample, in Feynman's book 19], where a minimal set of 
physical principles which must be satisfied in any theory 
of fields and particles is analyzed. Feynman includes in 
this set neither the second postulate of his formalism nor 
the Schrodinger equation: The only quantum mechanical 
principle included in this set is the principle of the su- 
perposition of probability amplitudes. In Ref. Q it has 
been shown that, instead of processes associated with a 
completely specified path in space-time, one can use pro- 
cesses associated with completely specified instants of the 
beginning and end of interaction in a quantum system as 
alternative ways in which any event can happen. As it 
turned out, employing this class of alternatives allows 
one to derive a dynamical equation from the principle 
of the superposition without making supplementary as- 
sumptions like the second postulate of Feynman's theory. 
By using this class of alternatives and the superposition 
principle, {ip\U{t,to)\f) ^ being the probability amplitude 
of finding a quantum system in the state {ip) in a mea- 
surement at time t, if at time to it was in the state \ip), 



can be represented in the form 

{^j\U{t,t„)\^) ^ {^\^) + f dt2 [ ' dt^{^\~S{hM)W).{l) 

Jto Jta 

Here (V'|>5'(i2, ii)|<y5) is the probability amplitude that, if 
at time ti the system was in the state \lp), then the in- 
teraction in the system will begin at time ti and end at 
time t2 and at this time the system will be in the state 
\'ip). By using the operator formalism, one can represent 
amplitudes {ili\U{t,tQ)\Lp) by the matrix elements of the 
unitary evolution operator U {t, to) in the interaction pic- 
ture. The operator S(t2,ti) represents the contribution 
to the evolution operator from the process in which the 
interaction in the system begins at time ti and ends at 
time ^2- As has been shown in Ref. Q, for the evolution 
operator Q to be unitary for any t and to the operator 
S{t2,ti) must satisfy the equation 

{t2-tl)S{t2,ti) 

= [ dti j dt3{t4-t3)S{t2,t4)S{t3,ti). (2) 
Jti Jti 

A remarkable feature of this equation is that it works as 
a recurrence relation and allows one to obtain the opera- 
tors S{t2jti) for any ti and t2, if S{t'2,t[) corresponding 
to infinitesimal duration times r = — t'l of interaction 
arc known. It is natural to assume that most of the con- 
tribution to the evolution operator in the limit t2 — > ti 
comes from the processes associated with the fundamen- 
tal interaction in the system under study. Denoting this 
contribution by Hint{t2,ti) we can write 

S{t2,ti) H,nt{t2,ti)+0{T'), (3) 

where t = t2 — ti. The parameter e is determined by 
demanding that Hint{t2,ti) called the generalized inter- 
action operator must be so close to the solution of Eq. 
(0) in the limit t2 — > ti that this equation has a unique 
solution having the behavior (|3J) near the point t2 = <i. If 
Hint{t2,ti) is specified, Eq. Q allows one to find the op- 
erator S'(t2,ti), and hence the evolution operator. Thus 
Eq. ^ which is a direct consequence of the principle 
of the superposition can be regarded as an equation of 
motion for states of a quantum system. This equation 
allows one to construct the evolution operator by using 
the contributions from fundamental processes as build- 
ing blocks. In the case of Hamiltonian dynamics the fun- 
damental interaction is instantaneous. The generalized 
interaction operator describing such an interaction is of 
the form 

H^nt{t2,ti) - -2i5{t2 - ti)Hi{ti) (4) 

(the delta function S{t2 —ti) emphasizes that the interac- 
tion is instantaneous). In this case Eq. ^ is equivalent 
to the Schrodinger equation Q and the operator Hj (t) is 
an interaction Hamiltonian (in the interaction picture). 
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At the same time, Eq. ^ permits the generahzation to 
the case where the fundamental interaction in a quan- 
tum system is nonlocal in time, and hence the dynamics 
is non-Hamiltonian 4]. 

By using Eq. for U{t,to), we can write 



U{t,ta) = 1 + — / dx eiip[-i{z - HQ)t] (5) 

X (z - Ho)-^T{z){z - Hoy' exp[i(z - Ho)to], 

where z = x + iy, y > 0, and Hq is the free Hamiltonian. 
The operator T{z) is defined by 



/•OO 

T{z) = i / dr exp{izT)T{T), 
Jo 



(6) 



where T(t) = exp(— ii/oi2)<5'(i2, ^i) exp(ii/oti), and t = 
t2 — ti. In terms of the T-matrix defined by Eq. jnj the 
equation of motion ^ can be rewritten in the form 0| 



d{ij2\T{z)\^i 
dz 



E 



{i^2\T{z)\n){n\T{z)\ij, 



(7) 



where n stands for the entire set of discrete and continu- 
ous variables that characterize the system in full, and |n) 
are the eigenvectors of H^- As it follows from Eq. Q), 
the boundary condition on this equation is of the form 

{^l^2\T{z)\^,) ^ {4>2\B{z)\i,{)+o{\z\-f') 

\z\ — > OO 

^ mB{z)\i^i) + O {h{z)} , (8) 

B{z)^i dTeM^zT)H\:l{T) (9) 
Jo 



where 



e^p{~iHot2)Hint{t2,ti) exp(ii7oii) 



is the interaction operator in the Schrodinger picture, 
and h{z) is an arbitrary function satisfying the condi- 
tion h{z) — o{\z\~^), \z\ OO, with l3 — e + 1. In 

the case of the Hamiltonian dynamics when H^^^{t) = 
^2i5{T)Hi, with Hi being the interaction Hamiltonian 
in the Schrodinger picture, B{z) — Hi, and the boundary 
condition IHl) takes the form 



T{z) 



Hi. 



(10) 



Equation Q with this boundary condition is equivalent 
to the Lippmann-Schwinger (LS) equation with the in- 
teraction Hamiltonian Hi. By definition, the operator 
B{z) represents the contribution which H-^^{t) gives to 
the operator T{z), and is the interaction operator in the 
energy representation. This operator must be so close 
to the relevant solution of Eq. |(7|) in the limit |2;| ^ cx3 
that this differential equation has a unique solution hav- 
ing the asymptotic behavior (|SJ). For this the operator 
B{z) must satisfy the condition 



d{iP2\B{z)\^i 
dz 



^ {i;2\B{z)\n){n\B{z)\^j,) 



(z-E, 



n) 



o[\z 



OO. 



(11) 



From Eq. lO it follows that the evolution operator in 
the Schrodinger picture can be represented in the form 



2ti 



dx exp{—izt)G{z), 



where z = x + iy, y > 0, and 



Giz) = Goiz) + Go{z)T{z)Go{z) 



(12) 



(13) 



with Go{z) — ^J^ii^ ■ Being equivalent to representation 
(0, Eqs. H12|l and (|13() express the principle of the su- 
perposition of probability amplitudes. In the ordinary 
quantum mechanics the similar equation plays an im- 
portant role and establishes the connection between the 
evolution operator and the Green operator G{z) which is 
defined by 



G{z) 



1 



H' 



(14) 



with H being the total Hamiltonian. Such a form of 
the Green operator follows from the fact that in the 
Hamiltonian formalism the evolution operator satisfies 
the Schrodinger equation. In the canonical formalism 
the T-matrix is defined by Eq. H13|l starting with the 
Green operator of the form (|14|l . In the formalism of the 
GQD the T-matrix plays a more fundamental role. It is 
defined by Eq. ©, and in this case the starting point is 
representation with S{t2,ti) being the contribution 
to the evolution operator from the processes in which 
the interaction begins at time ti and ends at time t2 ■ In 
this case the operator G{z) itself is defined by Eq. H13|l 
via the T-matrix. This is a more general definition of the 
Green operator, since representation (^) is a consequence 
of the principle of the superposition of probability ampli- 
tudes and must be valid in any case while the evolution 
operator can be represented in the form (|12(l with the 
operator G{z) given by Eq. (|14l) only in the case where 
the interaction in a system is instantaneous. 

As has been shown in Ref . , there is a one-to-one cor- 
respondence between the character of the dynamics and 
the large-momentum behavior of the T-matrix. If this 
behavior satisfies the requirements of ordinary quantum 
mechanics, then the interaction in a system is instan- 
taneous and the dynamics is Hamiltonian. In the case 
where this behavior is "bad", i.e., does not meet the re- 
quirements of the ordinary quantum mechanics, the in- 
teraction generating the dynamics of the system must 
necessarily be nonlocal-in-time. Let us now illustrate this 
point by using the model developed in Refs. 0, El as a 
test model demonstrating the possibility of the extension 
of quantum dynamics provided by the GQD. This model 
describes the evolution of the system of two identical non- 
relativistic particles whose interaction is separable, and 
hence the interaction operator has the form 



<P2|i/i:)(r)|pi >=^*(p2)^(pi)/(r). 



(15) 
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where p is the relative momentum of the particles, and 
/(r) is some function of the duration time r of the in- 
teraction in the system. It is assumed that in the limit 
IpI ^ oo the form factor </?(p) behave as 

'/'(p)-ipr", (IpHoo). (16) 

In this case the general solution of Eq. (U)) is 

'rf'*(P2)¥'(Pl) 



(p2|T(z)|pi) = — 



ga' + {z-a)Jj^j 



d^k |y(k)P 



,(17) 



z-Ek)(a-Ek) 



where ga — t{a), and a E (— oo,0]. In the case a > ^, 
the amplitude (p2|T(z)|pi) given by Eq. H17() . tends to 
A(y5*(p2)(^(pi), where 



Ek) ■ 

From Eqs. JHJ and © it follows that in this case the 
interaction operator H\'^^ (r) should be of the form 

< p2|ffi:)(T)|pi >= -2z,5(r)^*(p2)^(pi), 

and hence the dynamics of the system is Hamiltonian 
and is governed by the Schrodinger equation with the 
potential (p2|V'|pi) = \tf* {^2)f{T£>i) ■ 

In the case a < i, the T-matrix (|17|l tends to zero as 
\z\ oo 

(p2|T(z)|pi) ^ V*{V2MVl) 

\z\ ^ oo 

X (&i(-z)"-^ + 62(-z)2"-i) + o(|z|2"-i), (18) 

where bi = — 47r cos(a7r)m"^^ and 62 — 6i|a|^~" — 
bl{M{a)+g-^) with 



M{a) 



d^k |v?(k)p - |k|"^" 
(2^ a-Ek ■ 



In this case the dynamics is non-Hamiltonian and, as 
it follows from Eqs. © and ©, is generated by the 
nonlocal-in-time interaction operator 



<P2|lfi^)(T)|pi > = ¥'*(P2)^(Pl) 



X (aiT-"-^ -fa2r-2") , (19) 



where ai = Airi cos(a7r)m"~ 2 r~-^(i— a) exp[z(— -l-l-DTr], 
and 02 is a free parameter of the theory. The solution of 
Eq. with the interaction operator (|19|l is of the form 



(P2|r(z)|pi) = 



&?V'*(P2)V'(Pl) 



-62 + 6i(-z)5-"-M(z)6? 



(20) 



III. THE PIONLESS THEORY OF NUCLEAR 
FORCES 

The Weinberg program for low energy nucleon physics 
employs the analysis of time-ordered diagrams for the 2N 
T-matrix in ChPT to derive a NN potential and then to 
use it in the LS equation for constructing the full NN 
T-matrix. Obviously the starting point for this program 
is the assumption that in the nonrelativistic limit ChPT 
leads to low energy nucleon dynamics which is Hamilto- 
nian and is governed by the Schrodinger equation. How- 
ever, the fact that the chiral potentials constructed in 
this way are singular and lead to UV divergences means 
that this assumption has not corroborated. At the same 
time, the GQD allows one to analyze the predictions of 
ChPT without making a priori assumptions about the 
character of low energy nucleon dynamics: This charac- 
ter should results from the analysis. Let us consider, for 
example, the low energy predictions of ChPT for the 2N 
system in the ^5*0 channel. At very low energy, even the 
pion field can be integrated out, and the diagrams of the 
ChPT take the form of the diagrams being produced by 
the effective Lagrangian containing only contact interac- 
tions among nucleons and derivatives thereof |2Clj| . From 
the analysis of time-ordered diagrams of this theory it 
follows that the 2N T-matrix in the ^Sq channel must be 
of the form 



(p2|T(z)|pi) 



00 

E 



2n 2m, 

P2 Pi ''r. 



(21) 



where Pi is relative momentum of nucleons, and the 
terms tnm{z) are of order |t„™(z)| - 0{A-2(»+™)}, 
with the expansion parameter A being set by the pion 
mass. Obviously, chiral symmetry play no role in this 
case. However, for our analysis it is important that at 
extreme low energies ChPT gives rise to this result. 

From Eq. pi() it follows that at leading order the 2N 
T-matrix is of the form 



(p2|T(«)(z)|pi) = ioo(2). 



(22) 



On the other hand, this T-matrix must satisfy the gen- 
eralized equation of motion, and this requirement deter- 
mines its form up to one arbitrary parameter 



(P2|T(")(z)|pi) 



ra\J—zm 



(23) 



where m being nucleon mass. Note that the standard 
EFT approach yields the same expression for the lead- 
ing order T-matrix It is remarkable that the 
T-matrix (|23l) . which in the EFT approach is obtained 
by performing regularization and renormalization of the 
solution of the Schrodinger (LS) equation, is not a solu- 
tion of this equation. At the same time, it is a solution 
of Eq. (71) corresponding to the interaction operator |^ 

(P2l4:l WIPi) = 4«m-3/2 + Z7) , (24) 
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where 7 = 47rexp(— z7r/4)m~^/^C(^^. This operator is 
nonlocal in time and this is the only reason why in this 
case the generalized dynamical equation cannot be re- 
duced to the LS equation, and hence the T-matrix is not 
its solution. 

Let us now consider the full 2N T-matrix including all 
the terms in Eq. |j2lj. From Eqs. Q and it follows 
that the functions i,im(z) satisfy the equations 



dz 



(27r)3 ^ 

^ ^ n'm'=0 



(25) 

The solution of this equation is of the form 

tnyyiiz) = C2nC2mt[z) , (26) 

where t{z) satisfies the equation 



dt{z 



dz 



-t\z) 



d^k ?/>2(fc/A)V'i(fc/A) 
(2^ (z-i?fe)2 ■ 



(27) 



with 



oo 
n=l 

The form of the functions '0i(fc/A) and 4'2{k/A) mani- 
fests the fact that the constants C2n and C2„ are of order 
0{A^^"}, and hence the relevant domain of the defini- 
tion of these functions corresponds to the scale A. As has 
been shown in Ref. Q , for the T-matrix being a solution 
of Eq. H25(l to be unitary, it must satisfy the condition 

(P2|r+(z)|pi) = (p2|r(z)|pi), z e (-00, 0). 

From this and Eq. (|27|l it follows that 

Now we can rewrite Eq. (127(1 in the form 

di{z) f d^k 1 f d^k |V'(fc/A)|2-l 



dz J (27r)3 [z - Ekf J (2^)3 (z - Euf 
where t{z) — t~^{z) and 

oo 

V'(fc/A) = 7Ai(fc/A) = ^C2„fc2", 



(28) 



with Co = 1. The solution of this equation with the 
boundary condition <(z = 0) = to is 



i{z) — tQ — zm? 



d^k 



1 



(27r)3 (zm- fc2)fc2 



M(z), 



where 



M{z) = zm 



2 j d^k Fi{k/A) 
(27r)3 (zm - fc2)fc2' 



with Fi{k/A) = |V'(fc/A)|2 - 1. Hence 



t{z) 



47r 



\/-zm - M{z) 



where Co — . Thus from the requirement that the 2N 
T-matrix having the structure H21|l satisfy the generalized 
dynamical equation it follows that this T-matrix should 
be of the form 



(P2|T(z)|pi 



r{P2/A)^{pi/A) 
-'-^v^-M{zy 



(29) 



IV. THE EXPANSION OF THE EFFECTIVE 
THEORY. 

As we have shown, the requirement that low energy nu- 
cleon dynamics be governed by the generalized dynamical 
equation places additional constraints on the form of the 
2A^ T-matrix. This requirement together with the re- 
quirement that in the ^Sq channel the 27V T-matrix have 
the structure yield the expressions for this T- 
matrix which can be used for constructing the scattering 
amplitude and the evolution operator. The only prob- 
lem is that we do not know the values of the parameters 
C2n- One may hope that in future it will be possible to 
derive them from QCD. Another way is to obtain these 
parameters from low energy experiment. Of course, in 
this way one has to restrict oneself to a few parameters 
that should be fitted to experiment. In other words, one 
has to use Eq. H29|) in the spirit of the effective theory 
of nuclear forces, keeping only terms in the theory up to 
some order in Q/A in calculations of observables up to 
the corresponding order. 

In the EFT approach the scattering amplitude, which 
is given by the equation 



Aip) 



47r 



1 



TO p cot 6 ^ ip 



(30) 



(here and below we focus on the ^Sq partial wave), is de- 
scribed in terms of an expansion in powers of momentum 
over the characteristic short-distance scale A. It is well 
known that the quantity p cot S has a nice momentum 
expansion for p < < A known as the The Effective Range 
Expansion (ERE): 



pcotS 



1 1^ 

-a + 2^ 



2{n+l) 



(31) 



The scattering length a can be arbitrary large, while the 
other length scales ro, ri, ... are assumed to all have 
natural sizes, r„ ~ l/A(2"+i)_ Substituting Eq. (pTT|l in 
Eq. 1(30(1 we can expand the scattering amplitude about 
p = 0: 



, 47ra.^ 
Aip) = 1 - tap + aro/2 - 

TO 

+ o(pVa^)]. 



a2)p2 (32) 
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The radius of convergence of this expansion is about 
l/|a|, and hence in the case \a\ >> 1/A, which is of rel- 
evance to reahstic NN scattering, the expansion shown 
in Eq. H31|l breaks down for momenta far below A. In 
this case, one needs to expand the scattering amplitude 
in powers of p/ A while retaining ap to all orders: 



Aip) 



An 



-[1 



ma + tp 
^2 



ro/2 



iro/2y 



a + ip 



ip 



This is the Kaplan-Savage- Wise (KSW) expansion of the 
effective theory ,20j. It is the expansion that we wish to 



reproduce starting with the expression for the T-matrix 
given by Eq. if^ . Using Eq. for the scattering 

amplitude 

Aip) = -{p\Tiz = Ep + tO)\p), 

we can write 



A{p) 



mp/AW 



C^' + ^ip~M{Ep 



(33) 



In order to perform the momentum expansion of this am- 
plitude, let us represent M{z) as follows 



M{z) 



(fk Fi(fc/A) 2 
— - — - — = — 7Tn 

(27r)3 {zm ~ k^)k^ 

OO / 
71=1 \ 



(Pk Fi(fc/A)^^2^3 [ d^k Fi(fc/A) 



(2^) 

m 
in 



k^ 



i=0 



(27r)3 {zm-k^)k^ 

\ OO 



(34) 



n=l 



r 



with 



Mn{z) = -{zm)"J„ + ^\/^zm(zm)"^C2jC2(„_i), 



8=0 



where 



(35) 



Jn = m 



(fk F„{k/A) 

(27r)3 fc2(n+l) 



mA' 



.2„+i f d^k F^{k) 



Here fc = fc/A, 



(27r)3 k^(n+l) 



n—1 i 



(36) 



F^ik/A) = |^(fc/A)p -Y.Y.^2JC2i^-,)k'\ (37) 

By using Eq. H34() . we can represent the scattering am- 
plitude in the form 



A{p) = - 



ri=0 



2n 



1 + 37*^ E C^2„p2n 

Tl = 



(38) 



where C2n are determined by the recurrence relation 



C2n — Cq I C2iC2(„_i) — ^2(71- 



Equation H38|l allows one to expand the scattering am- 
plitude as 



A{p) = - 



Co 



C2P^ 



^iC2fip^ C,p^ 



Comparing with Eq. H32() , this expression relate the cou- 
plings C2n to the low energy scattering data a, r„: 

47ra 27r 2 

Go = , O2 = — a To, 

TO TO 

TO \ 4 2a 
By using the expression 



1 



47r 



(39) 



(40) 



n=0 n=0 

which relates couplings C2n(^) to C2n(0) = C2n as 



C2n{p-) 



n-1 

C271 + E C2i{p)C2{r. 

1=0 



Eq. H38|l can be rewritten in the form 



Aip) 



E C2n{^^)p^ 

n=0 



l + t(^P + M) E C2n{f^)p^ 
n=0 



(41) 



8 



This equation is exactly what Kaplan, Savage, and Wise 
have derived |23| in the pionless theory by using the PDS 
substraction scheme with the substraction point — 
— /z^, and Eq. H41|l yields the following expressions for 
Co (a*) and C2(/i): 



1 Tfl 



C2(m) 



C2 



Cn 



1/a — 

47r ro/2 
m (1/a — /i)^ 



As for the couplings C2n{fJ-) for n > 2, from Eq. H4Q|I it 
follows that they satisfy the equation 



^ m=0 



n 



i.e., satisfy the complete, coupled RG equation. 

We have demonstrated that the expression for the 2N 
T-matrix shown in Eq. H29(l can be used for calculating 
the scattering amplitude up to some order in Q/A. As 
we have seen, Eq. (|29ll reproduces the same expansion 
of the scattering amplitude that in the standard pionless 
EFT is derived by summing the bubble diagrams and 
performing regularization and renormalization. On the 
other hand, in general one cannot restrict oneself only 
to the scattering amplitudes. In the next section we will 
show that the formalism of the GQD provides a new way 
to formulate the effective theory of nuclear forces as a 
perfectly consistent theory that allows one to construct 
not only scattering amplitudes but also the off-shell 2A^ 
T-matrix, and the Green and evolution operators. 



V. THE EFFECTIVE THEORY OF NUCLEAR 
FORCES. 

As follows from Eq. in order to describe the low 
energy dynamics we have to obtain the T-matrix for any 
z relevant for the low energy theory. Equation (0) allows 
one to obtain the operator T[z) for any z provided that 
some boundary condition on this equation is specified. 
The boundary condition ^ means that the most of the 
contribution to the operator T(z) in the limit [zj ^ 00 
comes from the processes associated with the fundamen- 
tal interaction in the system under study that are de- 
scribed by the operator B{z). On the other hand, for the 
low energy theory to be consistent, the operator B{z) 
must be determined in terms of low energy degrees of 
freedom. This means that, despite the boundary condi- 
tion |(SJ formally implies that z must be let to infinity, 
really one has to restrict oneself to a "high" energy re- 
gion which is much above the low energy scale but is well 
below the scale of the underlying high energy physics. 
We will use V to denote this "high" energy region. The 
fact that the energy region V lies much above the scale of 
the low energy dynamics implies that at such " infinite" z 
the most of the contribution to the operator T(z) comes 



from processes that can be thought of as a "fundamen- 
tal" interaction in the low energy theory, and, as a result, 
the interaction operator B{z) is close enough to the true 
2N T-matrix. Thus, in order to take into account the 
fact that any theory has a range of validity, instead of 
the boundary condition ((HJ , we have to use the following 
boundary condition on Eq. ( |7| l: 



(V'2|T(z)|Vi 



= {^2\B{z)\i;{) 

+ 0{h{z)}, zeV. 



(42) 



For Eq. Q with the boundary condition l|42|l has a 
unique solution the interaction operator must be close 
enough inside domain V to the true 2N T-matrix. This 
means that it must satisfy the equation 



d(^2\B{z)\^i 
dz 



E 



{^2\B{z)\n){n\B{z)\^, 



[z - i?„)2 

+ O {\z\-^h{z)] , zeV. (43) 

In the low energy theory of nuclear forces all processes 
that in ChPT are described by irreducible diagrams in- 
volving only two external nucleons can be considered as 
such " fundamental" processes. Here irreducible diagrams 
are 2A^ irreducible: Any intermediate state contains at 
least one pion or isobar. It is natural to expect that the 
relative contribution of reducible 2N diagrams tends to 
zero as z increases, and, in a region that lies much above 
the scale of low energy nucleon dynamics, the main con- 
tribution comes from the "fundamental" processes that 
are described by the 2A^ irreducible diagrams, i.e., the op- 
erator B{z) becomes very close to the true 2A^ T-matrix. 
On the other hand, because of the separation of scales 
provided by QCD the above " high" energy region of the 
low energy theory lies still much below the scale of the un- 
derlying physics. In other words, the interaction operator 
and amplitudes {'4'2\B{z)\'il'i) that constitute the interac- 
tion operator and generate low energy nucleon dynamics 
really are low energy (in the scale of the underlying the- 
ory) amplitudes that in ChPT are described by the irre- 
ducible 2A^ diagrams. In principle they can be obtained 
within the underlying high energy theory in some low en- 
ergy limit. It is important that the matrix elements of 
the NN interaction operator are just the amplitudes that 
are described by these diagrams, and hence there is no 
a need to make a priory assumption that what we have 
to extract from the analysis of the diagrams in ChPT is 
the NN potential. The amplitudes that are described by 
these diagrams can be directly used as building blocks 
for constructing the low energy theory. Thus the GQD 
allows one to build a bridge between QCD and low en- 
ergy nucleon dynamics. It is hoped that in the future it 
will be possible to obtain the amplitudes {ip2\B{z)\il;i) in 
terms of QCD with such accuracy that the corresponding 
operator B(z) will determine a unique low energy 2N T- 
matrix. For this the operator B{z) obtained in this way 
must be close enough to the true T-matrix inside the 
domain V. There is no reason to believe that this opera- 
tor must necessarily generate the Hamiltonian low energy 



9 



dynamics. In fact, as it follows from Eq. H10() . for this 
(P2|-B(2:)|pi) must have a negligible dependence on z in- 
side the domain T). However, as we will show below, such 
a behavior of (p2|i3(z)|pi) is at variance with the sym- 
metries of QCD. On the other hand, QCD must produce 
the low energy 2iV T-matrix that satisfies the general- 
ized dynamical equation. From this in turn it follows 
that QCD must give rise to the interaction operator 
B{z) satisfying the condition 

In order to demonstrate how Eq. and the analysis 
of diagrams in ChPT can be used for constructing the 
interaction operator, let us turn to the pionless theory. 
As has been noted, from this analysis it follows that the 
2N T-matrix in the ^iSo channel must be of the form (|21|l . 
Obviously the interaction operator that reproduces 
this T-matrix inside the domain 2? must have the same 
structure 



(p2|B(z)|pi> - ^ 



2n 2m7 
Pi Pi ^nrn 



(z). 



(44) 



we have to keep the first two terms in the expansion of 
the function /(z) shown in Eq. H47I) . This means that the 
function h{z) in Eq. (|^ is determined by the equation 

h[z) ^0{\C,^M-\z)\] =0{\C,^M,\z)\] . (48) 
From the above it follows that the operator describing 
the NN interaction in the ^5*0 channel can be written as 

(p2|B(z)|pi> - -^*(p2/A)^A(Pi/A) (49) 
X {M-^{z) + Co^M-^{z)) . 

Correspondingly, for the interaction operator h\'^^{t), we 
have 

X (/i(r) + Co-V2(r)), 

where 



Correspondingly from the requirement that the operator 
B{z) have such a structure it follows that this operator 
should be of the form 



/»(!-) = — / dTexp(-izT)A^ "(z). 



(p2|B(z)|pi> = V/(P2/A)^(pi/A)/(z), 



(45) 



where ip{p/ A) is given by Eq. 
is a solution of the equation 



and the function /(z) 



dfjz) 
dz 



with 



J(z) = 



+ 0{\z\-^h{z)} , zeV, 
d'^k |V(fc/A)P 



(46) 



(27r)3 (z - Euf ■ 



Solving Eq. yields 

f{z) = -M-\z)-C^^M-\z) 
~Ca'^M-^{z)... + 0{h{z)}, zeP, 



(47) 



where M{z) ^ Mo(z) + M{z), with Mo(z) = ^^/^zm. 
Obviously the number of the terms we have to keep in 
this expansion depends on the value of the parameter /3, 
i.e., on the accuracy with which the operator B{z) of the 
form (|44|l reproduces the relevant solution of Eq. 10 in 
the limit |z| — > oo. Since the function A4{z) is completely 
determined by the form factor ?/'(p/A), all solutions of 
Eq. 146|l for given form factor have the same first term 
in the expansion 147(1 . and only the second term distin- 
guishes the different solutions (they are fixed by the free 
parameter Co). Thus in order that Eq. lO with the in- 
teraction operator of the form (|44|) has a unique solution, 



As we have shown, the knowledge of the form factor 
in Eq. H45() completely determines the form of the NN 
interaction operator which in turn determines the dy- 
namics of the system. Equation ((JJ) with this interaction 
operator has a unique solution, and this solution is the 
T-matrix shown in Eq. (|29|l . However, this implies that 
we know the details of the dynamics at high energies, 
because the relevant domain of definition of the func- 
tion ip{p/A) spreads to such energies. At the same time, 
the basic premise of effective theories is that low energy 
dynamics can be described in terms of a few parame- 
ters without any knowledge of the details of high energy 
interactions. As we have seen in Sec. IV, the functions 
ijj{p/A) and M{z) that enter the expression for the oper- 
ator B{z) can be expanded as is shown in Eqs. (|28|l and 
(|34|l . and in order to calculate these functions and hence 
(p2|-B(z)|pi) with accuracy up to the order (Q/A)^^ it is 
sufficient to know the first N terms in these expansions. 
In other words, for determining the interaction operator 
with the accuracy up to the order (Q/A)^^, one need 
only to know the properties of the form factor that are 
parameterized by the constants C2n, •■■ C2a^ and JJi, 
J^f^. Let us now show that, starting with the interac- 
tion operator (|49|l . where only the first N terms in the 
expansion of ^(p/A) and the same in the expansion of 
M(z) are known, one can calculate the T-matrix with 
the accuracy up to the order (Q/Ay^ for all relevant z. 
For calculating with this accuracy, instead of the inter- 
action operator (|49|l . one can use the following effective 
interaction operator 
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{p2\B^^]{z)\vi) = - (1 + c*^pI + ■■■ + c^pI-^ + v:kiP2/^)) (1 + + . . . + C2^pf^ + wiPi/^)) 



(50) 



with 



AT 



ri=0 
TV 



A/" 



— V-zm (zm)" ^ c;iC2(„_j) - ^(2:to)"J'„ 



n=0 i=0 



where if^(p/A) is an arbitrary function that satisfies the 
condition 

m{p/^) = o{{p/Af^^'^}. (51) 

The form of this effective interaction operator manifests 
uncertainty in the details of the NN interaction at high 



energies. In contrast with the operator B{z) which un- 
equaUy determines the relevant solution of Eq. 10, the 

effective interaction operator ^e// (•^) determines the set 
of the solutions that coincide inside the domain V with 
the true 2N T-matrix with the relative accuracy up to 
the order (g/A)2(-^+i) 



(P2|r(z)|pi) = (p2|i?(f;(z)|pi) (52) 
(l + 0(Q/A)2(^+i)) + 0{\Ca^M^'{z)\}, zeV. 



This implies that the operator b'^^J (z) must satisfy Eq. 
(UI) with the same accuracy 



J 



rf(P2ii3y;(z)ipi) 

dz 



d^k {p2\Bl^}{z)\k){k\Bi'^j.{z)\p^) 



(27r)3 (z - 

X (l + 0(Q/A)2(^+i)) + 0{\C^^\z\-^M,-^{z)\}, z e V. 



(53) 



By taking into account Eqs. (35-37), one can easily verify that this condition is really satisfied, provided the function 
(fij^ip/ A) is such that 



d^k 



|l + C2fc2 



C2^k^^ + W(fc/A) -EE 4,c2U-j)k^' 



i=Q j=Q 



(27r)3 



I — 



(54) 



where 1 < n < J\f. Thus the requirement that inside 
the domain T) the effective interaction operator is close 
enough to the true 2A'^ T-matrix determines some proper- 
ties of the function ipj^{p/A) that is the unknown part of 
the form factor A). The functions ^pj\f{p/ A) mainly 
characterize the form factor at high energies. For this 
reason, for describing low energy nucleon dynamics with 
the accuracy up to the order 0{Q / AY'^-^^^^ one need not 
to know the function ipj^{p/ A) exactly. This is because 
the details of the high energy physics must not affect on 
the low energy dynamics. It is sufficient to know the in- 
tegration properties of this function shown in Eq. (|54|l . 
Thus the constants J7i, J^j^ in the effective interaction 
operator B^j{z) parameterize the effects of high energy 
physics on low energy nucleon dynamics. 

By definition, the effective interaction operator 
B^^j{z) is not so close inside the domain V to the true 



27V T-matrix to determine a unique solution of Eq. Q. 
However, it determines a set of the solutions T'^^)(z) of 
Eq. lO that coincide with the true 2A'' T-matrix, with 
the relative accuracy up to the order 0(Q/A)2(^+i) (we 
will denote this set Each of these solutions corre- 

sponds to some definite function ipj\f{p/ A) satisfying the 
conditions H51() and (|54|l and can be represented as 

(p2|tW(z)|pi) = tMz)(Y.4y,- + ^*^{p2/A)] 



\n=0 



P?" + ^^'(Pi/A) , 



where t^{z) is the solution of the equation 



dz 



-tlf{z)J{z) = t^iz) 



dMjz) 
dz 



(55) 
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with the boundary condition 

t^iz) = - (Mj^i(z) + C^'Mj^^z)) (56) 
X (l + 0{{Q/Af^+'^}) + O {\C,H4,-Hm , 



z e V. 



The general sohition of Eq. (|55|l is 

1 



C-i-7W(z)' 



(57) 



where C is some arbitrary constant. Equation ()57|l may 
be rewritten in the form 



1 



C-^~MMz)-M^(z)^ 



where Mj^{z) = ^ Mn{z) is the part of Ai{z) that 

n=J\f 



depends on the choice of the function (pj^{p/A). From 
this equation it foUows that 



(l + 0{(Q/A)2(^+i)}) + O {\C-Hd^%z)\} , zeV. 



Here we have used the fact that Aij^{z) is of order 
O |((5/A)^''^^^'' |. On the other hand the relevant so- 
lution of Eq. H55|l must satisfy the boundary condition 
(|56|l . From this it follows that the parameter C in Eq. 

must be equal to Cq. Thus the solution of Eq. {T)) 
with the boundary condition given by Eqs. I|5UI) and 1521) 
where (pj\f(j)/A) is specified in some way is 



P2|r(^)(z)|pi) 



E cLpI" + '^^(P2/A)) ( E C2„p?" + w(Pi/A) 

n=0 / \n=0 

C^' - M^riz) - MMz) 
E cLpI" + ^*^{P2/A)) ( E C2„P?" + W(pi/A) 

n=Q / \n=0 

Co' - Mm{z) 

I 



(l + o{(Q/A)^(^+^)} 



(58) 



However, as has been noted, the function ipj\f{p/A) that 
enters the effective interaction operator is not known ex- 
actly. At a given order, only its integral properties shown 
in Eq. H54(l are known. For this reason Eq. H58|l actu- 
ally represents the set ^l_\f oi the solutions of Eq. ^ 
which coincide each with other with accuracy up to the 

order O |(Q/A)^^^+^^|, and the true 2N T-matrix be- 
longs to this set. This is because the "operator" B^^j{z) 
actually is a set of the interaction operators B{z) which 
coincide with the true 2N interaction operator with the 
same accuracy. This uncertainty in specifying the bound- 
ary condition on the generalized dynamical equation is 
only difference of the effective theory from a " fundamen- 
tal" one whose predictions are assumed to be exact. It 
is important that in both of these cases we deal with the 
same well defined equation of motion which can be used 
not only for analytic, but also for numerical calculations. 

As an example, let us consider the effective pionless 
theory at next-to-leading order {M =1). In this case, for 
the effective interaction operator, we have 



where 



(P2|i3S(^)|pi) 



= - 1 



* 2 
C2P2 



(1 



-^i(pi/A)) {M^\z) 



^\{P2/A)) (59) 
+ C^^M-,^{z)), 



Miiz) = — \J-zm - zmJi (-zm)^^'^Rec2, 

4:TT 2tt 

and ipi{p/A) is some function satisfying the condition 
shown in Eq. ifMll . The solution of Eq. ^ with this 
effective interaction operator yields 

(p2|T(i)(z)|pi) = (1 + c*pI + vlip2/A)) (60) 

{1+C2pl + Vl{pi/A)) 



Mi{z) 



[1 + 0{{Q/Ar}) 



This equation represents the set of the solutions of the 
generalized dynamical equation that coincide with the 
true 27V T-matrix with the relative accuracy up to the 
order (Q/A)^. At this order one may use any of these 

solutions. Correspondingly the operator BgJ:j(z) with 
any function ipi{p/A) satisfying the conditions 151(1 and 
(|54|l may be used as the interaction operator B{z) that 
determines a unique solution of Eq. Q. For exam- 
ple, one may chose this function as follows ipi{p/A) — 
C2P^exp(-p2/A^), with Ai = 2^T^/^^Jl{mRec2y^■ 

It should be noted that the expression (|58|l which rep- 
resents the form of the effective 2N T-matrix in the ^Sq 
channel of the pionless theory can also be derived from 
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Eq. (|29|l by keeping only the first Af terms in both the ex- 
pansions oftp{p/A) and Al{z) that enter in this equation. 
This is because in the case of the pionless theory the gen- 
eralized dynamical equation can be solved exactly, and as 
we have shown, the unique solution consistent with sym- 
metries of QCD is the T-matrix H29(l . However, this is not 
the case when pions must be included as explicit degrees 
of freedom. In the theory with pions, one cannot solve the 
problem exactly, and hence the way of constructing the 
27V T-matrix presented in Sec. Ill is inapplicable. Above 
we have presented the way of constructing the effective 
theory that prescribes to start with the analysis of the 
theory at "high" energies belonging to the domain T). 
At these energies the theory has the most simple struc- 
ture. For example, in the case where the dynamics in the 
theory is Hamiltonian, (-02 1^(^)1-01) have a negligible de- 
pendence on z inside the domain V. This means that 
the interaction generating the dynamics can be assumed 
to be instantaneous and hence the low energy dynamics 
can be described by an interaction Hamiltonian that is 
just the value of the T-matrix in this domain. If we deal 
with several different interactions, then the interaction 
Hamiltonian is the sum of the corresponding interaction 
Hamiltonians. This is because processes that involve sev- 
eral different interactions cannot occur instantaneously, 
and hence give a negligible contribution to the T-matrix 
inside the domain V. In the general case, however, the 
operator T{z) have significantly depend on z even inside 
this domain. As it follows from Eqs. I|29(l and (|58(l . this 
takes place in the case of the contact NN interaction. 

As an illustration of the advantages of the above for- 
mulation of the effective theory, let us consider the dy- 
namics of the 2N system in the presence of an external 
potential. For reason of physical transparency, we will 
consider nucleons as spinless particles whose interaction 
is described by the interaction operator H49() . i.e., is the 
same as in the ^Sq channel of the pionless theory. The 
interaction operator generating the dynamics of such a 
system should be of the form 



(P2;P2|BP°*(^)|Pi;pi) 



(P2;P2|S(z)|Pi;Pi) 
(P2;P2|V^|Pi;Pi), (61) 



where the vector |P; p) describes the state of the 2A^ 
system with the momentum of the center of mass of the 
system P and the relative momentum of nucleons p. The 



operator B{z) that describes the interaction of nucle- 
ons with themselves and, in the center of mass frame, 
is shown in Eq. 1)49(1 . The second term on the right hand 
side of Eq. (|61l) describes the instantaneous interaction 
with the external potential. We will assume that this 
potential is local, i.e., has the form 

(P2; P2|l^|Pi; Pi> = V(g)(27r)3^(3)(p^ - Pi - q) 

with q = P2 — Pi, and is week enough to be considered 
in the Born approximation. 

Let the details of the NN interaction are known only 
up to next-to-leading order. In this case we have to solve 
the problem in the spirit of the effective theory starting 
with the effective interaction operator of the form 



P2;P2|i?^;}(^)|Pi;Pi) 



(P2;P2|i?5^/(z)|Pi;pi) 
-l/(g),5(3)(p2-pi-q),(62) 



where ^^^■'^(z) is given by Eq. H59() . This operator deter- 
mines a set of the solutions of Eq. {Tj) that coincide with 
the true T-matrix describing the dynamics of the system 
with accuracy up to the order {Q/A)*. It can be checked 
that in the Born approximation these solutions are given 
by the equation 

TP°*(z) =T(i'(z) 
+ (1 + T(i)(z)Go(z)) ^ (1 + Go(z)T(i)(z)) , (63) 

where T(i)(z) is one of the solutions of Eq. shown 
in Eq. H60() . By using Eq. H63() . one can construct, 
with a given accuracy, the scattering amplitude and the 
evolution operator. For example, for the amplitude of 
2N scattering in the presence of the external potential, 
we can write 

A(P2,P2;Pi,Pi) = -(P2;P2| (1 + T(i)(z)Go(z)) (64) 
xy(l + Go(z)T«(z)) |Pi;pi), 



where z = -r — h — -l-iO = -r^ + — + iO, and pi ^ 
P2- By using Eq. ((60(l . for the scattering amplitude 
A(P2,p2;Pi,Pi) we get 



J 



A(P2,P2;Pi,Pi) = yloo(P2,P2;Pi,Pi) +^oi(P2,P2;Pi,Pi) + ^io(P2,P2;Pi,Pi) + ^ii(P2,P2;Pi,Pi), (65) 
where 

Aoo(P2, P2; Pi, Pi) - -F(g)(27r)35(3)(p2 - p^ - q), 



A..(P.P.P..P.) = - '^""^l^":"-^-'l-' (l + 0(Q/An 



Zi-Ep- 



Pi 
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V{q){l + clp^ + c^pl){l + 0{Q/KY) 
y(g)(l + c^pi+C2p'^)(l + 0(Q/A)4) 



An(P2,P2;Pi,Pi) = -^(9) 



7 (27r)3 (z-iJp, -£;p^+jO)(z-£;p-i;p, +zO) ^ ^ J 

^ ng)(l+C^Pi+C2pf) 



(27r)3 (i;p, - i;^/ + iQ){Ep, - Ep + iO) J (27r)3 Ep{Ep., - i;^, + iO)(£:p, - Ep + iO) 



Here we have used the fact that the function 931 (p/A) sat- 
isfies the condition (|54|) . Thus the generahzed dynamical 
equation with the effective NN interaction operator H50() 
ahows one to calculate, with the desired accuracy, not 
only 2N data but also processes in the presence of an ex- 
ternal field. As we have seen, in this way regularization 
and renormalization are not required, and hence one need 
not to introduce new parameters: The solution depends 
only on the constants Co, C2, and J7i that are contained in 
the effective interaction operator. As we have noted, the 
coupling Co is fixed by the scattering length Co — ^a. 
However, knowing the 2A^ scattering data is not sufficient 
to obtain the parameter C2 and J7i. Only the coupling 
C2 = Co(c2 + c* - Co Ji) = Co2§a is fixed by the scatter- 
ing data at this order, but not the couplings C2 and J7i 
separately. At the same time, as it follows from Eq. (|66|l . 
these parameters manifest themselves in 27V scattering in 
the presence of an external potential. This is because the 
2N T-matrices enter half-of-the-energy shell in the scat- 
tering amplitude Ht)4|) describing such processes. In other 
words, such processes can be used for obtaining the above 
couplings. 



VI. THE PROBABILISTIC FRAME OF 
QUANTUM MECHANICS AND THE 
CHARACTER OF LOW ENERGY NUCLEON 
DYNAMICS. 

In an EFT renormalization allows one to describe low 
energy physics using an effective Lagrangian that con- 
tains only a few degrees of freedom, because after renor- 
malization that consists in absorbing infinities in a re- 
definition of constants in the Lagrangian the remaining 
integrals are effectively cut off at internal momenta lager 



than Q. This might lead us to the conclusion that the 
definition of any effective theory must necessarily include 
regularization and renormalization. However, as we have 
shown, by using the example of the pionless theory, reg- 
ularization and renormalization are needed when we de- 
scribe the theory in terms of an effective Lagrangian. 
But, if we consider the problem from the more general 
point of view provided by the GQD and do not restrict 
ourselves to the assumption that the effective action is 
instantaneous, we see that the effective theory of nuclear 
forces manifests itself as a perfectly consistent theory free 
from the UV divergences. As we have shown, in the pio- 
nless theory the 2A^ T-matrix in the ^5*0 channel is given 
by Eq. H29(l . This T-matrix satisfies the generalized dy- 
namical equation, i.e., the solution of this equation does 
not require regularization and renormalization. In other 
words, in the pionless theory the generalized dynamical 
equation allows one to separate the low energy physics 
from the underlying high energy physics without renor- 
malization. Let us now show that this directly follows 
from the first principles of quantum mechanics. 

The fact that because of the separation of scales, at low 
energies nucleons emerge as the only effective degrees of 
freedom means that the temporal evolution of the system 
can be described in terms of the 2N Hilbert subspace 
Ti.2N and the evolution operator defined on this space 

U2N{t,0)^P2NUs{t,0)P2N, (67) 

where P2N being the projection operator on the subspace 
H2N- Here we use the Schrodinger picture. Since in 
this case the 2N system is considered to be closed, the 
evolution operator C/2Ar(t, 0) should be unitary 

U+j^{t,0)U2Nit,0) =1. (68) 
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This is one of the main requirements that quantum 
mechanics imposes on such a theory. It expresses the 
fact that, if at initial time the system was in the state 
IV') G 'H2N, then in a measurement at time t the system 
wiU be necessarily found in one of the 2N states. Here, 
of course, it is assumed that H2N describes 2N states at 
low energies. The above means that the sum of the prob- 
abilities to find the system in all the possible 2N states 
must be equal to unite 



(2^ 



\{k\u2N{t,om\' ^1, 



(69) 



with lijj) being a normalized vector belonging to TL2N- 
Here we consider nucleons as spinless particles, and 
|(k|C/2Ar(<2,*i)IV')P(|:;;^ is the probability of finding the 
quantum system in the 2N state with the relative mo- 
mentum p in the volume of the momentum space ki < 
Pz < ki + dk, [i = 1,2,3). Finally, for Eq. ^ to be 
valid for any normalized vector \^) g Ti.2N, the operator 
Us{t,0) must be unitary. 

Another basic principle of quantum mechanics is the 
principle of the superposition of the probability ampli- 
tudes from which it follows that the evolution operator 
is defined by Eq. H12|) which in this case reads 

(p2|t/2Ar(i,0)|pi) = ^ J dxexp{-izt)G2N{z), (70) 



with 



where 



G 



2n{z) = G'.iz) + G'o{z)T2n{z)G',{z), (71) 



G'oiz) ^ P2nGo{z)P2n 
and T2n{z) is defined by the equation 



(72) 



00 

(P2|72w(2)|pi) = i J d{t2 - h) exp[iz{t2 - h)] 



X (p2|T2Ar(t2-il)|pi), 



(73) 



Thus, in the case of the reduced 2N system, we deal 
with the same equations as in describing the complete 
system. The only difference is that the amplitudes 
(P2|'Z2Ar(<2 — ^i)|Pi) describes the contributions to the 
evolution operator from the processes in which the inter- 
action in the 27V system (not in the complete system) 
begins at time ti and ends at time ^2- By using Eq. (|70|l . 
the unitarity condition (|68|l may be rewritten as 



oc oc 



(27r 



^ J dxi J dx2exp[it{z2 - zi)] 



-00 —00 



For Eq. (|7Hl to be valid for any time t, the 27V T-matrix 
must satisfy Eq. iQ. Note that in Ref. the general- 
ized dynamical equation has been derived just in this way. 
At the same time, the above does not mean that T2n{z) 
satisfies the LS equation. In other words, from the ex- 
istence of the separation of scales which results in the 
unitarity of the 2A^ evolution operator U2N{t, 0) and the 
superposition principle it follows that the 2iV T-matrix 
must necessarily satisfy the generalized dynamical equa- 
tion but not the LS equation. Thus the dynamics of the 
2iV system may be non-Hamiltonian. 

Let us now show that the above is true even if the 
dynamics in the underlying theory is governed by the 
Schrodinger equation. As we have noted, in the case 
where the dynamics of a quantum system is Hamiltonian 
the evolution operator Us{t, 0) can be represented in the 
form H12() with the Green operator given by Eq. H14|l . 
From this it follows that 

U2N{t,0) ^ P2NUs{t,0)P2N 

i f°° 

= — / dxexp{-izt)P2Niz - H)-^P2 

Stt 7_oc 

/oo 
dxexp{-izt)G2N{z)- (75) 
-oc 

It can be shown, that the Green operator G2N — 
P2n{z — H)~^P2N may be represented in the form 



2N 



G2n{z)^P2n{z-H)-^P2n 
= G',{z) + G'^{z)T2n{z)G'^{z), 



(76) 



where Go(z) is defined by Eq. H72II . The above means 
that, if the underlying dynamics is governed by the 
Schrodinger equation, then the 27V evolution operator 
U2Nit,0) can be represented in the form 



t/2Af(i,0) ~ — / dx exp{—izt) 
X iG'oiz) + G'o{z)T2N{z)G'„{z)). 



(77) 



x{'4j2\G+^{z2)G2n{ziMi) = {H^i)(7A) 



From this it follows that T2n{z) must satisfy Eq. (0, be- 
cause the operator (|77|l is assumed to be unitary. At the 
same time, the above does not mean that the amplitudes 
(k2|T(z)|ki) fall off rapidly enough above the scale of the 
effective theory for the dynamics in the theory to be gov- 
erned by the Schrodinger equation. In fact, the Hamilto- 
nian character of the underlying dynamics implies only 
that the 2A^ T-matrix must fall off rapidly above the 
scale of the underlying physics. Thus, the dynamics in 
the effective low energy theory may be non-Hamiltonian 
even if the underlying high energy dynamics is governed 
by the Schrodinger equation. 

An important lesson to be learned from the above anal- 
ysis is that low energy nucleon dynamics can be described 
by the generalized dynamical equation reduced to the 
2iV Hilbert space without regularization and renormal- 
ization. The divergence problems may arise only in de- 
scribing the dynamics of the effective theory by using 
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the Schrodinger (LS) equation, because such a descrip- 
tion is based on the additional assumption that the ef- 
fective NN interaction is instantaneous and, as a conse- 
quence, the generahzed dynamical equation is equivalent 
to the Schrodinger equation. Note, in this connection, 
that from the physical point of view the fact that the 
chiral potentials lead to UV divergences means that the 
LS equation with such potentials does not provide a well 
separation of the low energy physics from the underly- 
ing high energy physics. In other words, the amplitudes 
(k2|r(2;)|ki) does not fall off rapidly enough above the 
scale of the effective theory for integrals in the LS equa- 
tion to be effectively cut off at internal momenta larger 
than Q. In this case, in order to integrate out the high 
energy degrees of freedom, one need to use regularization 
and renormalization. In contrast, the generalized dynam- 
ical equation should separate the low energy physics that 
is described by the effective theory of nuclear forces from 
the underlying high energy physics provided that at very 
low energies nucleons may emerge as the only explicit de- 
grees of freedom. In fact, the letter means that the matrix 
elements of the operator U2N{t,0) must fall off rapidly 
enough at momenta larger than Q. This is because the 
probability of finding, for a measurement at time t, the 
system in the 2N state with momenta larger than Q, if 
at time t = the system was in the low energy state 
must be negligible. If this is not the case, one can- 
not also ignore other states with such momenta in which 
the high energy degrees of freedom manifest themselves, 
and hence the unitarity condition (|69|) breaks down. As 
it follows from Eqs. 0, and the fact that 

(k2|J7(<, to)|ki) fall off rapidly at momenta larger than 
Q in turn means that the integrals in the generalized dy- 
namical equation effectively cut off at such momenta, and 
hence, in this case the high energy degrees of freedom can 
be integrated out without renormalization. 



VII. A NEW LOOK AT THE WEINBERG 
PROGRAM 

Our formulation of the effective theory of nuclear forces 
may be regarded as a new way of realizing the Weinberg 
program for physics of the two-nucleon systems. The im- 
plementation of this program has three stages. Firstly, 
one must consider the case where the NN system is not 
subject to any external probes and employ ChPT to gen- 
erate a nonrelativistic particle-number-conserving Hamil- 
tonian for the nuclear system. At the second stage this 
Hamiltonian should be employed for constructing the full 
2N T-matrix via the LS (the Schrodinger) equation. The 
third step of the Weinberg program is to employ the 2N 
T-matrix obtained in this way in the case where 2N sys- 
tem is subject some external probe: a pion, a photon, 
or some weakly-interacting particle. Provided this probe 
carries momentum of order m^r its interaction with the 
2N system can be represented as a sum of irreducible 
diagrams which forms a kernel Kprobe for the process of 



interest. The full amplitude for this process is then found 
by multiplying this kernel Kprobe by the factors describ- 
ing the interaction of the NN pair in the initial and final 

states 

A={1+ TGo)Kprobe{l + GoT). (78) 

The Weinberg proposal was based on the assumption 
that the only equation that can govern low energy nu- 
cleon dynamics is the LS (the Schrodinger) equation, and 
hence what one has to derive from the analysis of dia- 
grams in ChPT is an NN potential. However, there is 
no reason to consider that low energy nucleon dynamics 
is necessarily governed by the Schrodinger equation. In 
principle this dynamics may be governed by the general- 
ized dynamical equation with a nonlocal-in-time interac- 
tion operator when this equation is not equivalent to the 
Schrodinger equation. In fact, as has been shown in Ref. 
4] , only the generalized dynamical equation must be sat- 
isfied in any case, not the Schrodinger equation. In the 
light of this fact the Weinberg program can be considered 
from a new point of view: Instead of the Schrodinger (LS) 
equation, one should use the generalized dynamical equa- 
tion. In this case one need not to use a priory assumption 
that from the analysis of diagrams for the 2N T-matrix 
it follows that the an effective NN interaction is instan- 
taneous. Such a modification of the Weinberg program 
is natural and does not change its character. Indeed, 
the fact that the analysis of the diagrams in ChPT leads 
to singular chiral potentials in the case of which the LS 
equation makes no sense without renormalization means 
that really this equation is not sufficient for construct- 
ing the 2N T-matrix. In addition one need to specify 
a renormalization scheme to make these predictions fi- 
nite. But the LS equation plus a renormalization scheme 
is something more general than this equation itself, and, 
after renormalization we deal with the dynamics that is 
governed by another equation which, as has been shown 
in Ref. Q , may be only the generalized dynamical equa- 
tion. In order to illustrate this point let us consider the 
effective theory at leading order of the Weinberg power 
counting. The NN effective potential that has been de- 
rived from the Weinberg analysis of diagrams in ChPT 
is ^1] 

+ C5 + CtCi • 172, 

with q = P2 — Pi. The coupling is the axial cou- 
pling constant, is the pion mass, /jr is the pion de- 
cay constant, and (t(t) are the Pauli matrices acting in 
spin (isospin) space. At very low energies the one-pion- 
exchange part of the NN interaction may be included 
into the contact term, and the NN potential takes the 
form 

(p2|y|pi) = C, (79) 
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where Cq = Cs — 3Ct + g\/ f^. Obviously, the potential 
(|79|l is singular, and the LS equation with this poten- 
tial makes no sense without regularization and renor- 
malization. Let us perform regularization by using a 
momentum cut-off. In this case the singular potential 
(|7(^ is replaced by the regularized one (P2|Va|pi) = 
/*(p2/A)Co(A)/(pi/A),where the form factor f{p/A) 
satisfies /(O) = 1 and falls off rapidly for p/A > 1. The 
solution of the LS equation with this potential is 



P2|rA(z)|pi) =r(p2/A)/(pi/A) 
d^k \f{k/A)\^'-' 



X C,\A) 



(2^)3 z-Ek 

Defining the renormalized value Cr of Co as the value of 
the T-matrix (pa |7a(z)|pi) atz = ^ = ^= 



d'k \f{k/AW 



Ek 



(27r)3 

we may rewrite this equation as 

/*(p2/A)/(pi/A) 



(80) 



(P2|rA(^)|Pi) 



c; 



d^k |/(fc/A)P 
{z-E^)Ek 



where the renormalized value is fixed by the scattering 
length a = mCi^/47r. Thus, after renormalization the 
integral in the expression for the leading order 27V T- 
matrix is effectively cut off, and hence at this stage the 
regularization may be removed by letting A — > cx), and 
we get 



(P2|r(z)|pi)= c^^(A)- 



^3/2 

47r 



(81) 



Correspondingly, for the half-of-the-energy shell T ma- 
trix, we have 



(p'|T(i?, + iO)|p)= 



-1 , *"^lPl 



47r 



This is just what has been obtained by Weinberg in Ref. 
P . Obviously the T-matrix is not a solution of the 
LS equation. At the same time, as we have seen, this T- 
matrix is the well defined solution of the generalized dy- 
namical equation with the interaction operator H24II de- 
scribing a nonlocal in time interaction. The dependence 
of this interaction operator on z expresses the fact that 
the instantaneous interaction which is described by the 
contact potential (|79|) does not contribute to the 2N T- 
matrix separately. In fact, as it follows from Eq. H80() af- 
ter removing regularization, this potential becomes equal 
to zero. A nonzero contribution to this T-matrix comes 
only from the sum of diagrams involving infinite number 
of such contact interactions. In other words renormaliza- 
tion spreads effective contact NN interaction in time. 

Thus even if we start with the LS equation and try to 
extract an NN potential from the analysis of diagrams 



in ChPT, despite it is singular and makes no sense, after 
renormalization we come to low energy nucleon dynamics 
that is governed by the generalized dynamical equation 
with a nonlocal- in-time interaction operator. This means 
that in order to realize the Weinberg program in a con- 
sistent way one has to use this equation instead of the LS 
equation. As we have shown in the previous section, in 
contrast with the LS equation, the generalized dynamical 
equation separates the low energy physics from the high 
energy physics without renormalization. This manifests 
itself in the fact that, being reduced to describing the 
27V system, this equation is free from UV divergences. 
In Sec.V we have demonstrated the advantages of this 
way of realizing the Weinberg program by using the ex- 
ample of the pionless theory. Firstly, by using the analy- 
sis of diagrams in ChPT, we have obtained the effective 
7V7V interaction operator. In the channel it is of the 
form 150|) . By solving the generalized dynamical equa- 
tion with this interaction operator we have constructed 
the ^Sq channel 27V T-matrix shown in Eq. It is 

important that in this way we have obtained not only the 
27V scattering amplitude but also the off-shell T-matrix. 
Equation (|58|l expresses the constrains that the symme- 
tries of QCD place on the off-shell behavior of the 27V 
T-matrix. At a given order there are only a few free pa- 
rameters in Eq. (|58|l that can be derived from low energy 
experiment. Once they are obtained in any way, the off- 
shell 27V T-matrix that in this case is completely deter- 
mined by Eq. (|58|l can be used in Eq. (|54|) for describing 
the interaction of the 27V system with an external probe. 
In Sec.V we have demonstrated this fact by using the ex- 
ample of the 27V scattering in the presence of an external 
potential. As we have seen, with the 27V T-matrix H58|) 
in hand the full amplitudes of the processes involving the 
interaction with external probes can be calculated with- 
out regularization and renormalization. 

The key point of the above way of the realization of 
the Weinberg program is the derivation of the 7V7V in- 
teraction operator from the analysis of the time-ordered 
diagrams in the time-ordered diagrams in the energy re- 
gion v. This region is sufficiently above the scale of the 
low energy physics for the processes that are described by 
the 27V T-matrix at 2; g P may be considered as a " funda- 
mental" interaction in this theory. At the same time, this 
energy region is much below the scale of the underlying 
high energy physics, and hence the effective 7V7V inter- 
action operator is really determined by the low energy 
in the scale of QCD behavior of the time-ordered dia- 
gram for 27V T-matrix. Moreover, one need not to know 
this behavior exactly. It is sufhcient to know the struc- 
ture of this T-matrix, i.e., its dependence on momenta 
of nucleons, that is predicted by ChPT. For example, in 
the pionless theory the channel 27V T-matrix has the 
structure shown in Eq. (|21|l . Correspondingly the struc- 
ture of the interaction operator is given by Eq. H44() . On 
the other hand, this interaction must be close enough in 
the energy region T) to the dynamical equation, i.e., must 
satisfy Eq. H43|) . This requirement yields the expression 
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for the effective interaction operator shown in Eq. (|45|l . 
This interaction operator can then be used for construct- 
ing the 2N T-matrix. At the same time, in the pionless 
theory the requirement that the 2N T-matrix having the 
structure H21() satisfy the gencrahzed dynamical equation 
directly yields the expression for the 2N T-matrix shown 
in Eq. (|29|l . and this expression can be used for organiz- 
ing the calculations of the effective 2N T-matrix. This 
is because the 2N T-matrix has the same simple struc- 
ture at all relevant z as in the "high" energy region T). 
However, this is not the case in the theory with pions 
or (and) when one must take into account the Coulomb 
interaction between protons. In such theories the struc- 
ture of the 27V T-matrix is much simpler in the "high" 
energy region than at low energies where the nonpertur- 
bative character of nucleon dynamics manifests itself. In 
this case only in the region V the 2N T-matrix, with the 
accuracy needed for constructing the interaction opera- 
tor, can be represented as a sum of contributions from 
the contact, pion-exchange, and Coulomb interactions. 
The fact that starting from this high energy representa- 
tion one can construct the NN interaction operator and 
then use them for obtaining the 2N T-matrix has been 
demonstrated in Ref. (2^] by using the example of the 2N 
dynamics at leading order of the Weinberg power count- 
ing. As we have seen, the same situation takes place in 
the case when one consider the 2A^ dynamics in presence 
of an external potential. In this case interaction opera- 
tor generating the dynamics in the system is given by Eq. 
()61|l . and its form manifests the fact that the "high" en- 
ergy region V the main contribution to the 2N T-matrix 
can be represented as a sum of the contributions from the 
interaction of nucleons with themselves and the external 
potential. By using the interaction operator H62|l one can 
describe the dynamics of this system. 



VIII. THE OFF-SHELL BEHAVIOR OF THE 2iV 
T-MATRIX AND THE THREE-NUCLEON 
PROBLEM 

As we have shown, the formulation of the effective the- 
ory of nuclear forces presented in Sec. V allows one to 
construct not only the 27V scattering amplitude (in this 
case we reproduce all results of the standard EFT ap- 
proach), but also the off-shell T-matrix, and the evolu- 
tion and Green operators. This is very important be- 
cause the S'-matrix is not everything. For example, at 
finite temperature there is no S'-matrix because parti- 
cles cannot get out to infinite distances from a collision 
without bumping into things. This means that the off- 
shell structure of the 2iV T-matrix should influence on in- 
mcdium observables. In Ref. j8j it has been shown that 
the off-shell behavior of the in-medium nucleon-nucleon 
T-matrix and hence the off-shell properties of nuclear 
forces have substantial effects on the transition ampli- 
tudes and cross sections at large nuclear matter densities. 
These properties are crucial for solving the many-nucleon 



problem. For example, the 2iV amplitudes enter off-the- 
energy-shell in the 3A^ equations. 

The realistic NN potentials that describe 27V scatter- 
ing data to high precision can not guarantee that a simi- 
lar precision will be achieved in the description of larger 
nuclear systems. In fact, the simplest observable in the 
3A^ system, the binding energy of the triton, is under 
predicted by the realistic NN potentials which are so 
successful in describing the 2A^ observables. The energy 
deficit ranges from 0.5 to 0.9 MeV and depends on the 
off-shell and short-range parameterization of the 2iV force 
[l^ . In order to resolve this problem one has to take into 
account three-nucleon force (3NF) contributions to the 
ZN binding energy. The common way of solving the 3A^ 
bound state problem is to use in the Schrodinger equation 
phenomenological NN potentials and then to introduce 
a 3NF to provide supplementary binding. However, from 
the point of view of the three-nucleon problem, it is not 
sufficient to generate a phenomenological NN potential 
that perfectly reproduces the 2A^ scattering amplitudes. 
One must also generate a NN potential by using theo- 
retical insight as much as possible in order to constrain 
the off-shell properties of the 2A^ T-matrix. If this is 
not the NN potential which fits precisely the 2A^ 

phase shifts but produces the erroneous off-shell behav- 
ior of the T-matrix would not provide reliable results for 
the 3A^ system, nor can be used to test for the pres- 
ence of 3A'^ forces. It is important, in this context, that 
the requirement that the 2N T-matrix satisfy the gen- 
eralized dynamical equation and has the form consistent 
with the symmetries of QCD places a constrain on the 
off-shell behavior of this T-matrix. As we have shown, 
in the pionless theory the 2 T-matrix in the ^ Sq chan- 
nel should be of the form l|29() . The expression shown in 
Eq. (|29|) contains the infinite set of parameters C2„ that 
in principle could be obtained in terms of QCD. At the 
same time, as we have seen, the generalized dynamical 
equation allows one to organize calculations of the two- 
nucleon T-matrix in the spirit of the effective theory by 
parameterizing the effects of the underlying physics in a 
few parameters that can be derived from a low energy 
experiment. In this way we have obtained the expression 
(|58|l that determines 2 A*" T-matrix in the ^S'o channel 
with the accuracy up to the order (Q/A)^'^^^^ . For ex- 
ample, at next-to-leading order the T-matrix is given by 
Eq. H6UI) containing only three free parameters J7i , Rec2 
and lmc2 (the parameter Co is fixed by the scattering 
length at leading order). One of them is fixed by the 2N 
phase shifts analysis while other two can be derived from 
the in-medium scattering data or from the low energy 
data corresponding to the processes of the interaction of 
the two-nucleon system with an external probe. In other 
words, the requirement that the 2A^ T-matrix be con- 
sistent with basic principles of quantum mechanics and 
the symmetries of QCD removes the off-shell ambigui- 
ties. This may provide a better understanding of the 
many-nucleon problem. 

The models for the "iNY that are usually used for solv- 
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ing the problem with the triton understanding are based 
on two-pion exchange with intermediate A-isobar excita- 
tion. However, these 3A^F models cannot explain the Ay 
puzzle. In Ref. |3| a three-nucleon force generated by 
the exchange of one pion in the presence of a 2N corre- 
lation ^183 has been suggested as a possible candidate to 
explain the Ay puzzle. This 3NF contribution is fixed 
by the 2N T-matrix describing the underlying 2N inter- 
action while the pion is " in flight" . The expression for 
this force derived in Ref. jl^ contains the off-shell 2N 
T-matrix, more precisely its subtracted part 

il2(P2,Pl,2;) = tl2(P2,Pl,2;) - Wl2(P2,Pl), (82) 

where pi and p2 are Jacobi momenta of nucleons 1 and 
2, and the potential - like term Vi2(p2, Pi) contains only 
OPE/OBE - type diagrams. The subtraction in Eq. 
is needed to take into account a cancellation effect which 
has been observed in Refs. |2^ |2^. This cancellation 
involves meson retardation effects of the iterated Born 
term, and the irreducible diagrams generated by sub- 
summing all time ordering diagrams describing the com- 
bined exchange of two mesons amongst the three nucle- 
ons. In principle there are no free parameters to adjust, 
and the 3A^ force is completely determined by the 2N 
T-matrix. However, as has been shown in Ref. |2^, in 
order to explain the Ay puzzle, instead of the subtracted 
T-matrix shown in Eq. H82|l . one has to use the amplitude 
defined according to the prescription 

il2(P2,Pl,-Z) = c(z)ti2(p2,Pl,2) - Wl2(p2,Pl), 

with the effective parameter c(z), which represents an 
overall correction factor for the far-off-the-energy-shell 
2N T-matrix. Ideally, this parameter should be one for 
the 2N potential to provide a reliable extrapolation of 
the 2N T-matrix down to z « —160 MeV. However, as 
has been shown in Ref. none of the existing 2N 

T-matrices can guarantee the off-shell behavior that is 
needed for the explanation of the Ay puzzle with the 
parameter c(z) set to one. For example, in order to re- 
produce the nd experimental data with the Bonn B po- 
tential, the factor c{z) must be set to 0.73 for the energy 
3 MeV. 

The lesson one must learn from the above is that the 
ofF-shell behavior of the 2N T-matrix may play a cru- 
cial role in explaining the Ay puzzle, and the existing 
realistic potentials do not provide the off-shell behavior 
that is needed for the correct reproduction for Ay with 
the 3A^F suggested in Ref. Ts'l. Note in this connection 
that, being fitted to the 2N scattering data, the exist- 
ing realistic potentials lead to ambiguities in the off-shell 
behavior of the 2N T-matrix. In order to remove the off- 
shell ambiguities one has to find the way of constructing 
the 2N T-matrix as an inevitable consequence of the ba- 
sic principles of quantum mechanics and the symmetries 
of QCD. The remarkable feature of the EFT approach is 
that it implies to build the theory of nuclear forces as a 
consequence of these principles. However, the standard 



EFT of nuclear forces does not predict the off-shell be- 
havior of the 2N T-matrix. The cause of this is that 
in this theory the Schrodinger (LS) equation makes no 
sense without regularization and renormalization. On 
the other hand, there is no reason to restrict ourselves 
to the assumption that low energy nucleon dynamics is 
governed by the Schrodinger equation: Only the gener- 
alized dynamical equation must be satisfied in any case, 
and, as we have seen, the requirement that the solutions 
of this equation that are consistent with the symmetries 
of QCD correspond to a nonlocal-in-time interaction op- 
erator when the generalized dynamical equation cannot 
be reduced to the Schrodinger equation. At the same 
time, as has been shown in Ref. the nonlocality in 
time of the interaction results in an anomalous off-shell 
behavior of the T-matrix. In other words, the formalism 
of the GQD predicts that in order that the 2N T-matrix 
to be consistent with the symmetries of QCD it must 
have the off-shell behavior that cannot take place in the 
case of the ordinary potentials. From this point of view, 
the " anomalous" off-shell behavior that the 2 T-matrix 
must have for the 3NF be able to explain the Ay puzzle 
may be considered as a manifestation of the fact that low 
energy nucleon dynamics is really non-Hamiltonian. 



IX. SUMMARY AND DISCUSSION 

We have shown that from the Weinberg analysis of 
time-ordered diagrams for the 2A^ T-matrix in ChPT it 
follows that nucleon dynamics at low energies is governed 
by the generalized dynamical equation with a nonlocal- 
in-time interaction operator. A remarkable feature of the 
generalized dynamical equation which follows straightfor- 
wardly from the first principles of quantum mechanics is 
that it allows one to construct all physical amplitudes 
relevant for the theory under consideration by using the 
amplitudes describing processes in which the duration 
time of interaction is infinitesimal. It is natural to as- 
sume that the most of contribution to these amplitudes 
comes from the processes associated with a fundamen- 
tal interaction in a quantum system. This point mani- 
fests itself in the boundary condition ||2J). If we do not 
consider a theory that is valid up to infinitely high ener- 
gies (infinitesimal times), then the above infinitesimal of 
the duration times of interaction should mean that be- 
ing much smaller than the scale of the theory these times 
may be much larger than the time scale of the underlying 
high energy physics. This in turn means that the ampli- 
tudes describing the " fundamental" interaction in the low 
energy theory can be computed in terms of the underly- 
ing high energy physics. Thus the generalized dynamical 
equation allows one to take into account that every the- 
ory with which we deal is a low energy approximation to 
a more fundamental one and provides a bridge between 
them. 

The generalized dynamical equation can be repre- 
sented in the form of the differential equation Q for the 
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operator T(z). The boundary condition on the general- 
ized dynamical equation of this form is shown in Eq. ^ , 
where the operator B{z) describes the fundamental inter- 
action in this system. By definition, this operator must 
be so close to the true T- matrix in the limit |z| ^ cx) that 
the generalized dynamical equation with the boundary 
condition have a unique solution. The above means 
that really this region of "infinite" energies with which 
we have to deal is the domain 2? that lies much above 
the scale of the low energy dynamics but much below the 
scale of the underlying high energy physics. Correspond- 
ingly for describing the low energy dynamics we have to 
start with the boundary condition H42|l that implies that 
the interaction operator B{z) is so close to the relevant 
T-matrix inside the domain D that Eq. with this 
initial condition has a unique solution. In the theory of 
nuclear forces this domain is a region of energies that 
are high enough for the most of contribution to the 2N 
T-matrix to come from processes that are described by 
the irreducible 2N diagrams for the 2N T-matrix (these 
processes are associated with a "fundamental" interac- 
tion), but not so high for the heavy degrees of freedom 
manifest themselves explicitly. By using the analysis of 
the time-ordered diagrams for the 2N T-matrix in ChPT 
inside the domain D where the structure of the theory is 
much simpler than in the low energy region, we can ob- 
tain the NN interaction operator which should be used in 
the boundary condition (|42|l on Eq. (|7|) . We have shown 
that the dynamics which is generated by the interaction 
operator obtained in this way is non-Hamiltonian. This 
is because, for low energy nucleon dynamics to be Hamil- 
tonian, the operator T(z) must have a negligible depen- 
dence on z inside the domain V, while, as we have shown, 
this is not the case. In other words, in the nonrelativistic 
limit QCD leads through ChPT to low energy nucleon 
dynamics that is not governed by the Schrodinger equa- 
tion. However, this does not mean that the low energy 
predictions of QCD are not consistent with quantum me- 
chanics. This means only that in this case we deal with 
a nonlocal-in-time interaction when the generalized dy- 
namical equation cannot be reduced to the Schrodinger 
equation. 

We have shown that from the fact that at extreme low 
energies nucleons emerge as the only effective degrees of 
freedom it follows that the evolution operator U2Nit,to) 
defined on the 2N subspace 'H2N is unitary and can be 
represented in the form H73() . This in turn means that 
the 27V T-matrix that enters Eq. H74(l must satisfy the 
generalized dynamical equation. In other words the dy- 
namics in a consistent effective theory of nuclear forces 
must be governed by the generalized dynamical equation, 
and the problems with UV divergences can arise only if 
we restrict ourselves to the assumption that the effective 
NN interaction is instantaneous and hence the general- 
ized dynamical equation is equivalent to the Schrodinger 
equation. What we have to do then is to obtain the inter- 
action operator B{z), i.e., the operator that is so close to 
the true solution that the generalized dynamical equation 



with the boundary condition (jSJ has a unique solution. 
Ideally, this operator should be derived from QCD. How- 
ever, for describing low energy nucleon dynamics with a 
given accuracy one need not to know the operator B(z) 
exactly. It is sufficient to know the effective interaction 
operator B^^j (z) that determines the solution of Eq. (0 
with the same accuracy. The advantage of using such 
an interaction operator whose form manifests the sym- 
metries of QCD is that in this case the effects of high 
energy physics on low energy nucleon dynamics are pa- 
rameterized by a few constants. The effective interaction 
operator B^^j (z) determines the set ft^/ of the solutions 
of Eq. ITJl that coincide with the true two-nucleon T- 
matrix with the accuracy up to the order (Q/A)^"^, i.e., 
it represents the set of the interaction operators B(z) 
that determine the above solutions. This uncertainty in 
specifying the initial condition for the generalized dy- 
namical equation manifests the fact that we do not know 
the details of the underlying physics. The effects of this 
physics on the low energy dynamics are parameterized 
by the constants C2, . . ., C2^f and Ji, . . ., contain- 
ing in the operator Bf^j{z) that is shown in Eq. H5U|I . 
As N increases the set of the solutions of Eq. {T)) 
that are determined by the effective interaction operator 
B^^}{z) become smaller and smaller. This means that 
in this case one approaches to the true solution closer 
and closer, because the true 27V T-matrix enters each 
of these sets. Each of the T-matrices belonging to the 
set r^Tv" corresponds to some function ipj\f{p/A), which 
must be such that the effective interaction operator sat- 
isfies Eq. (|43|l . This requirement puts the constraints 
on the function (pf^{p/A): it must satisfy the condition 
(|54|l where the parameters jTn enter the effective inter- 
action operator through the function M„(z). In other 
words, the effective interaction operator contains all pa- 
rameters that appear in the theory at a given order, and 
these parameters need not be redefined in the process of 
calculations. From this point of view the only difference 
of the effective theory from a full one is uncertainty in 
specifying the boundary condition on the generalized dy- 
namical equation, and the boundary condition l|42l) with 
the effective interaction operator of the form (|50|l means 
that one may choose any function (pj\f{p/ A) satisfying the 
conditions H51|l and 1)54(1 . This is equivalent to the choice 
of the interaction operator B(z) that formally determines 
a unique solution of Eq. (jT)). Thus being formulated in 
terms of the GQD the effective theory of nuclear forces 
can be put on the same firm theoretical grounds as the 
quantum mechanics of atomic phenomena. In this case 
we deal with a well defined equation of motion that allows 
one to construct not only the scattering amplitudes but 
also the off-shell T-matrix, and the evolution and Green 
operators. 

Our formalism may be regarded as a new way of realiz- 
ing the Weinberg program for physics of the two-nucleon 
systems based on the use of the generalized dynamical 
equation instead the LS equation. In this case at the 
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first stage of the program we derive from the analysis of 
the diagrams ChPT an effective interaction operator. In 
the ^5*0 channel of the pionless theory this operator is of 
the form 150|) . Solving the generalized dynamical equa- 
tion with this interaction operator yields the expression 
for the 2N T- matrix shown in Eq. ISHJ- This T-matrix 
leads to the KSW expansion [2fl|] of the scattering ampli- 
tude shown in Eq. (|40|) . However, in contrast with stan- 
dard approach to the EFT of nuclear forces where this 
expansion is obtained by summing the bubble diagrams 
and performing regularization and renormalization, we 
reproduce the same result in a consistent way free from 
UV divergences. At the same time, Eq. (I55II determines 
the off-shell 2N T-matrix in the ^S'o channel up to a few 
constants that parameterize the effects of the underlying 
physics on low energy nucleon dynamics. In other words, 
the requirement that low energy nucleon dynamics be 
consistent with the symmetries of QCD leads to a con- 
strain on the off-shell behavior of the 2N T-matrix, and 
one of the lessons that one must learn from the results of 
our work is that at very low energies (in the pionless the- 
ory) the ^5*0 ofT-shell 2N T-matrix consistent with the ba- 
sic principles of quantum mechanics and the symmetries 
of QCD must, be of the form lf55|l . With this T-matrix in 
hand implementation of the third step of the Weinberg 
program can proceed without resorting to regularization 
and renormalization. For example, formula (|65l) we have 
derived allows one to calculate, in a consistent way free 
from UV divergences, the 2N scattering amplitudes in 
the presence of an external potential to subleading order. 
This potential, for example, may be the Coulomb poten- 
tial. Of course, in this case an additional term describ- 
ing the contribution from gauge-invariant, four-nucleon, 
one-photon operator that arise at next-to-leading order 



must be included. However, in our case this term does 
not play a role of a regulator, like the counterterms that 
are used in the standard EFT approach for calculating 
such observables as the electric-quadrupolc moment of 
the deuteron. It should be emphasized that, if the ex- 
pression shown in Eq. H66(l where divergent, then this 
term could not absorb divergences for all momenta of the 
nucleons. Moreover, the same off-shell 2N T-matrix can 
be used in describing the interaction of the 2N system 
with other probes, in the three-body continuum calcu- 
lations, and in solving the microscopic nuclear structure 
problems. The advantage of the formulation of the effec- 
tive theory of nuclear forces based on the GQD is that 
it allows one to construct the 2N interaction operator 
and hence the off-shell 2N T-matrix as inevitable con- 
sequence of the basic principles of quantum mechanics 
and the symmetries of QCD. This opens new possibili- 
ties for solving many problems in nuclear physics, and, 
in particular, for explaining the Ay puzzle. 

In this paper we focused on the pionless theory be- 
cause in this simple case the proposed way to formulate 
the effective theory of nuclear forces can be investigated 
in detail, and its ideas can be verified exactly. At the 
same time, the approach should be applicable to the the- 
ory with pions. The only problem is that in this case 
the definition of a consistent power counting scheme and 
the derivation of the effective operator of the NN inter- 
action from the analysis of the diagrams in ChPT are 
more complicated. However, once this operator is con- 
structed the generalized dynamical equation can be used 
for performing not only analytic but also numerical cal- 
culations of any observables without regularization and 
renormalization. 
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